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1. Introduction 

Quantum information aims at the treatment and transport of information using 
the laws of quantum physics. For these goals, continuous variables (CV) of the 
electromagnetic field have emerged as a powerful tool ^ |21 El- In this context, 
entanglement is an essential ressource. The purpose of this paper is to link the very 
powerful mathematical description of gaussian states based on covariance matrices 
and the experimental production and manipulation of entanglement. 

Experimentally, entanglement can be obtained directly by type-II parametric 
interaction deamplifying either the vacuum fluctuations as was demonstrated in 
the seminal experiment by Ou et al. @] (or in recent experiments 13 E)) or the 
fluctuations of a weak injected beam [7]. It can also be obtained indirectly by mixing 
on a beam splitter two independent squeezed beams. The required squeezing can 
be produced by Kerr effects - using optical fibers [H| or cold atoms in an optical 
cavity [H] - or by type-I parametric interaction in a cavity ^3 Single-pass 
type-I interaction in a non-colinear configuration can also generate directly entangled 
beams as demonstrated recently by Wenger et al. in the pulsed regime All 
these methods produce a symmetric entangled state enabling dense coding, the 
teleportation of coherent ^3 El ^3 or squeezed states ^1] or entanglement swapping 

I14| . These experiments generate an entangled state with a covariance matrix 
in the so-called 'standard form' ^| ^], without having to apply any local linear 
unitary transformations such as beam-splitting or phase-shifts to exploit it optimally 
in quantum information protocols. 

However, it has been recently shown in |2()j that, when a bircfringcnt plate is 
inserted inside the cavity of a type-II optical parametric oscillator, i.e. when mode 
coupling is added, the generated two-mode state remains symmetric but entanglement 
is not observed on orthogonal quadratures: the state produced is not in the standard 
form. The entanglement of the two emitted modes in this configuration is not optimal: 
it is indeed possible by passive non-local operations to select modes that are more 
entangled. Our original system provides thus a good insight into the quantification 
and manipulation of the entanglement resources of two-mode Gaussian states. In 
particular, as just anticipated, it allows to conflrm experimentally the theoretical 
predictions on the entangling capacity of passive optical elements and on the selection 
of the optimally entangled bosonic modes [21] . 

The paper is organized as follows: we start by giving a general overview of 
gaussian states and defining the covariance matrix formalism, focusing in particular 
on two-mode states. Over the last years, a great deal of attention has been 
devoted to defining not only criterion of entanglement but also quantification of this 
entanglement. Sectional focus on different such measures which are interpreted in the 
covariance matrix formalism. Effects of mode coupling on the generated entanglement 
is then discussed in detail (section 0)). The experimental setup used to generate 
and quantify the entanglement is presented in section El In section El experimental 
measures of entanglement are realized and we discuss on the effect of noise of the 
covariance matrix on the entanglement measures. Finally, we show how optimization 
of the resource is obtained by a passive operation - namely a polarization adjustment 
- operated on the two- mode state (section CJ. 
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2. Gaussian states: general overview 

A continuous variable (CV) system is described by a Hilbert space Ti = T^i 
resulting from the tensor product of infinite dimensional Fock spaces Tii's. Let 
and aj be the annihilation and creation operators acting on Tii (ladder operators), 
and Xi = {ui + a]) and pi = {ai — a\)/i be the related quadrature phase operators. 
The corresponding phase space variables will be denoted by Xi and pi. Let X = 
{xi,pi, . . . ,XmPn) denote the vector of the operators Xi and pi. The canonical 
commutation relations for the Xi can be expressed in terms of the symplectic form VL 

[Xi,Xj\ = 2i^lij , 

n 

with = , Lo = 

The state of a CV system can be equivalently described by a positive trace-class 
operator (the density matrix g) or by quasi-probability distributions. In the following, 
we shall adopt the Wigner quasi-probability function W{R) defined, for any density 
matrix, as the Fourier transform of the symmetrically ordered characteristic function 

W{R) = \ ! Tr [qDy] e'^^^^ d'^^Y , i? e R^n ^ 

where the displacement operators Dy (describing the effect on the field of a classical 
driving current) are defined as 

The 2n-dimensional space of definition of the Wigner function, associated to the 
quadratic form f2, is referred to as phase space, in analogy with classical Hamiltonian 
dynamics. In Wigner phase space picture, the tensor product H = of the Hilbert 

spaces Tii's of the n modes results in the direct sum F = of the phase spaces 

r,'s. 

States with Gaussian characteristic functions and quasi-probability distributions 
are referred to as Gaussian states. Such states are at the heart of information 
processing in CV systems ^ and are the subject of our analysis. By definition, 
a Gaussian state g is completely characterised by the first and second statistical 
moments of the field operators, which will be denoted, respectively, by the vector 
of first moments X = ^(Xi), (X2), . . . , (X2„-i), (X2„)^ and the covariancc matrix 
(CM) F of elements 

F,, EE i {X,X, + X,X,) - (X,) {X,) , (2) 

where, for any observable 6, the expectation value (6) = Tr{gd). Notice that, 
according to our definition of the quadrature operators in terms of the ladder operators, 
the entries of the CM are real numbers. Being the variances and covariances of 
quantum operators, such entries are obtained by noise variance and noise correlation 
measurements. They can be expressed as energies by multiplying them by the quantity 
hio, where lo is the frequency of the considered mode. In fact, for any n-mode state 
the quantity hwTr (F/4) is just the contribution of the second moments to the average 
of the "free" Hamiltonian Y^^=ii^i^i + 1/2)- 
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Coherent states, resulting from the apphcation of displacement operators Dx 
to the vacuum state |0), are Gaussian states with CM F = 1 and first statistical 
moments X = X. First moments can be arbitrarily adjusted by local unitary 
operations (displacements), which cannot affect any property related to entropy or 
entanglement. In the present experimental case, the fields are produced by an OPO 
below the oscillation threshold, and the means are actually zero. Let us note that this 
can be done without any loss of generality. 

The canonical commutation relations and the positivity of the density matrix g 
imply 

T + iQ>0, (3) 

meaning that all the eigenvalues of the (hermitian) matrix T + ifl have to be greater 
or equal than zero. Inequality Q is the necessary and sufficient constraint the matrix 
F has to fulfill to be a CM corresponding to a physical Gaussian state IS] ■ Note 
that the previous condition is necessary for the CM of any (generally non Gaussian) 
state, as it generalises to many modes the Robertson-Schrodinger uncertainty relation, 
reducing to the familiar Heisenberg principle for pure, uncorrelated states. We mention 
that such a constraint implies F > 0. 

A major role in the theoretical and experimental manipulation of Gaussian states 
is played by unitary operations which preserve the Gaussian character of the states 
on which they act. Such operations are all those generated by terms of the first and 
second order in the field operators. So, beside the already discussed displacements, 
also the unitary operations generated by terms of the second order are Gaussian. As 
a consequence of the Stone- Von Neumann theorem, any such operation at the Hilbert 
space level corresponds, in phase space, to a symplectic transformation, i.e. to a 
linear transformation S which preserves the symplectic form fi, so that ft = S'^flS, 
i.e. it preserves the commutators between the different operators. Symplectic 
transformations on a 2ri-dimensional phase space form the (real) symplectic group, 
denoted by Sp^2n,'SR.)- Such transformations act linearly on first moments and "by 
congruence" on covariance matrices {i.e. so that F i— > S'^TS). One has Det S* = 1, 
VS* G 5p(2„3). 

Ideal beam splitters, phase shifters and squeezers are described by symplectic 
transformations. In fact single and two-mode squeezings occurring, respectively, in 
degenerate and non degenerate parametric down conversions, are described by the 
operators 

f^»j,r,v = e^^™-"!-"*"'"^) with e = re'2'^,re]R,(/9e [0,2^], (4) 

resulting in single-mode squeezing of mode i for i = j. The representation in phase 
space of the operation Uij,r.ip for i 7^ j is given by the linear transformation Sij,r^ip 



/ c~hs ks \ 

c + hs -ks 

ks c + hs 

\ -ks c-hs J 



for i^ j , (5) 



where c = cosh(2r), s = sinh(2r), h ~ cos(2iy9), k ~ sm(2ip) (r and ip are same as in 
Eq. I0J) and the matrix is understood to act on the couple of modes i and j. Beam 
splitters are described by the operators 

Oy,e= e""^--'-^^ [0,2^], (6) 
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forz^j. (7) 



corresponding to symplcctic rotations Rij^ in phase space 

/ cos(6') -sin(6l) \ 

cos(6') -sin(6l) 

sin(6l) cos(6l) 

\ sin(6i) cos(6i) / 

The angle 6 is defined by Eq. © and. again, the matrix is understood to act on the i 
and j modes. Single-mode symplcctic operations are easily retrieved as well, being just 
combinations of single mode (two dimensional) rotations and of single mode squeezings 
of the form diag(e'', e^*") for r > 0. Now, symplectic transformations in phase space 
are generated by exponentiation of matrices written as Jf2, where J is antisymmetric 
|25|. Such generators can be symmetric or antisymmetric. The operations Rij,e, 
generated by antisymmetric operators are orthogonal and, acting by congruence on 
the CM r, preserve the value of TrF. Since TrF gives the contribution of the second 
moments to the average of the Hamiltonian aja^, these transformations are said to 
be 'passive', or 'energy preserving' (they belong to the compact subgroup of Sp(^2n,'^))- 
Instead, operations Sij^r,ip, generated by symmetric operators, are not orthogonal and 
do not preserve TrF (they belong to the non compact subgroup of Sp(^2n,'R))- This 
mathematical difference between squeezers and phase space rotations accounts, in 
a quite elegant way, for the difference between 'active' (i.e. energy consuming) and 
'passive' (i.e. energy preserving) optical transformations. 

Finally, let us recall that, due to Williamson theorem the CM of a n-mode 
Gaussian state can always be written as 

F = S^iyS , (8) 

where S € S'p(2ri.R) ^-^d v is the CM 

V = diag(z^i, 1^1, ... , i/„, i^n) , (9) 

corresponding to a tensor product of thermal states with diagonal density matrix £»® 
given by 

^^ = (8)^E(^)l^M^K (10) 

i * fc=0 ^ ^ 

\k)i being the fc-th number state of the Fock space Hi. The dual (Hilbert space) 
formulation of Eq. © then reads: g = W U, for some unitary U. 
The quantities Vi 's form the symplectic spectrum of the CM F and can be computed 
as the eigenvalues of the matrix \i^ir\ Such eigenvalues are in fact invariant under 
the action of symplcctic transformations on the matrix F. 

The symplectic eigenvalues Vi encode essential informations on the Gaussian state F 
and provide powerful, simple ways to express its fundamental properties. For instance, 
in terms of the symplectic eigenvalues Vi, the uncertainty relation ^ simply reads 

(11) 

Also the entropic quantities of Gaussian states can be expressed in terms of their 
symplectic eigenvalues and invariants. Notably, the purity Tr g^ of a Gaussian state g 
is simply given by the symplectic invariant Det F = Y[i=i being Tr g^ = (Det F)~^/^ 
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2.1. Two-mode states 

Since this work is focused on two-mode Gaussian states, we briefly review here some 
of their basic properties. The expression of the two-mode CM F in terms of the three 
2x2 matrices a, f3, 7 will be useful 



a 


7 


1' 


fi 



(12) 

For any two-mode CM F there exist local symplectic operations 5*1 and 5*2 (each Sj 
acting on one of the two modes), such that their direct sum Si = Si®S2 (corresponding 
to the tensor product of local unitary operations) brings F to the so called standard 
form Tsf mini 



Si TSi 



a 





c+ 


^ 





a 





c_ 


c+ 





b 








c_ 





b 



\ 



I 



States whose standard form fulfills a ^ b are said to be symmetric, 
recall that any pure state is symmetric and fulfills c+ = — c_ 



1 



(13) 



Let us 
The 



correlations a, b, c-|_, and c_ are determined by the four local symplectic invariants 



(ab — c\_){ab — c-_), Deta 



Det/3 = 6 , Det7 = c-|_c_. Therefore, the 



DctF 

standard form corresponding to any CM is unique (up to a common sign flip in c_ 
and c+). 

For two-mode states, the uncertainty principle Ineq. (jSjl can be recast as a 
constraint on the 5'p(4_r) invariants DetF and A(F) = Deta + Det/3 + 2Det7: 

A(F) < 1 + DctF . (14) 

The symplectic eigenvalues of a two-mode Gaussian state will be named v- and 
v^, with < i^+i with the Hcisenberg uncertainty relation reducing to 

A simple expression for the !/zp can be found in terms of the two 5'p(4_iB 
(invariants under global, two-mode symplectic operations) |29l I3L)| 



(15) 
invariants 



2vi 



,^ - A(F) T v/A(F)2 - 4 Dct F . (16) 

A subclass of Gaussian states with a major interest in experimental quantum 
optics and in the practical realization of CV quantum information protocols is 
constituted by the nonsymmetric two-mode squeezed thermal states. Let S^ = 
Si2,r,Tr/4 bc the two modc squeezing operator with real squeezing parameter r, and let 
£1®; be a tensor product of thermal states with CM v^,^ = l2^+, where is 
the symplectic spectrum of the state. A nonsymmetric two-mode squeezed thermal 
state ^y^^r is defined as ^^.,r = SrgfiS^, corresponding to a standard form with 



cosh r + sinh 



sinh 



. cosh r 



(17) 



± 



1^+ 



sinh 2r 



In the symmetric instance (with 1/^ = ly^ = v) these states reduce to two-mode 
squeezed thermal states. The covariance matrices of these states are symmetric 
standard forms with 



i^cosh2r, c± = ±i/sinh2r . 



(18) 
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These are the states which occur in most reahstic parametric down conversion 
processes, hke the one which wiU be discussed later on in the paper. In the pure case, 
for which ^ = 1, one recovers the two-mode squeezed vacua. Such states encompass 
all the standard forms associated to pure states: any two-mode Gaussian state can 
reduced to a squeezed vacuum by means of unitary local operations. 



3. Entanglement of Gaussian states 

This section aims at reviewing the main results on the qualification and quantification 
of entanglement for Gaussian states of CV systems, which will be exploited in the 
following. 

The positivity of the partially transposed state (PPT criterion) is necessary and 
sufficient for the separability of two-mode Gaussian states ^Hl (and, more generally, of 
all (1 +7T.)-mode Gaussian states under 1 x n,-mode bipartitions and of bisymmetric 
(to + ri,)-mode Gaussian states under to x n-mode bipartitions |32|). In general, 
the partial transposition ^ of a bipartite quantum state g is defined as the result 
of the transposition performed on only one of the two subsystems in some given basis. 
It can be promptly inferred from the definition of the Wigner function W{X) that 
the action of partial transposition amounts, in phase space, to a mirror reflection 
of one of the four canonical variables. In terms of S'p(2.R) © SP{2,R) invariants, 
this reduces to a sign flip in Det7. Therefore the invariant A(r) is changed into 
A(r) = A(f) = Deta + Dct/3 - 2Det7. Now, the symplectic eigenvalues of f 
read 



\ 



A(r) T A/A(r)2 -4Dctr 

^ \ • (19) 

The PPT criterion thus reduces to a simple inequality that must be satisfied by the 
smallest symplectic eigenvalue j>_ of the partially transposed state 

£>_>!, (20) 

which is equivalent to 

A(r) < Detr + 1 . (21) 

The above inequalities imply Det7 = c+c_ < as a necessary condition for a two- 
mode Gaussian state to be entangled. Therefore, the quantity v- encodes all the 
qualitative characterization of the entanglement for arbitrary (pure or mixed) two- 
modes Gaussian states. 

Let us now briefly focus on the entanglement qualification of symmetric states, 
which will be the subject of the experimental investigations presented in the paper. 
It is immediately apparent that, because a = b, the partially transposed CM in 
standard form F (obtained by flipping the sign of c_ ) is diagonalized by the orthogonal, 
symplectic transformation i?i2.ir/4 of Eq- 05 resulting in a diagonal matrix with 
entries a=F |cip|. The symplectic eigenvalues of such a matrix are then easily retrieved 
by applying local squeezings. In particular, the smallest eigenvalue v- is simply given 

by 



u. = V{a-\c+\)ia-\c^\). (22) 

Note that also the original standard CM P with a = b could be diagonalized {not 
symplectically, since the four diagonal entries are generally all different) by the same 
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beam splitter transformation i?i2,7r/4, with the same orthogonal eigenvalues |czp | . It 
is immediate to verify that is just given by the geometric average between the two 
smallest of such orthogonal eigenvalues of F. The two quadratures resulting from the 
previous beam splitter transformation select orthogonal directions in phase space with 
respect to the original ones, so they will be referred to as 'orthogonal' quadratures. 
Notice that, in the experimental practice, this allows for the determination of the 
entanglement through the measurement of diagonal entries (noise variances) of the CM 
after the application of a balanced beam splitter (which embodies the transformation 

-Rl2,7r/4)- 

To explore further consequences of this fact, let us briefly recall some 
theoretical results on the generation of entanglement under passive (energy-preserving) 
transformations, which will be precious in the following. As shown in Ref. )21j . 
the minimum value for i>_ {i.e. the maximal entanglement) attainable by passive 
transformations is given by ~ A1A2, where Ai and A2 are the two smallest 
eigenvalues of F. Therefore, the entanglement of symmetric states in standard form 
cannot be increased through energy preserving operations, like beam splitter and phase 
shifters. On the other hand, as it will be carefully discussed in great detail in the next 
section, the insertion of a birefringent plate in a type-II optical parametric oscillator 
results in states symmetric but not in standard form. In such a case the entanglement 
can be optimized by the action of a (passive) phase shifter. 

A measure of entanglement which can be computed for general Gaussian states 
is provided by the negativity A/", first introduced in Refs. PHll^j later thoroughly 
discussed and extended in Ref. to CV systems. The negativity of a quantum state 
Q is defined as 

M{,) = , (23) 

where q is the partially transposed density matrix and ||o||i = Tr|6| stands for the 
trace norm of the hermitian operator 6. The quantity M{q) is equal to | Ai|, the 
modulus of the sum of the negative eigenvalues of p, quantifying the extent to which 
Q fails to be positive. Strictly related to M is the logarithmic negativity £'^, defined 
as Ej^ = log2 ll^lli, which constitutes an upper bound to the distillable entanglement 
of the quantum state g. Both the negativity and the logarithmic negativity have 
been proven to possess the crucial property of being monotone under LOCC (local 
operations and classical communications) [221 12^] ■ 

For any two-mode Gaussian state g the negativity is a simple decreasing function 
of which is thus itself an inverse quantifier of entanglement: 



0: 



(24) 



Ej^{g) = max [0, - loga J>-] . (25) 

This expression quantifies the amount by which PPT inequality H20() is violated. 
The symplectic eigenvalue i>- thus completely qualifies and quantifies the quantum 
entanglement of a Gaussian state F. 

For symmetric Gaussian states one can also compute |19| the entanglement of 
formation Ep EHl- We recall that the entanglement of formation Ep of a quantum 
state g is defined as 

Epig) = min Y^pMA)) , (26) 

iPiAi'i)} , 
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where the minimum is taken over all the pure states realizations of g: 

Q = ^Pi\'^i){'^i\ ■ 

i 

The quantity Ep satisfies all the requirements of a proper entanglement measure. 
The asymptotic regularization of the entanglement of formation is equal to the 
entanglement cost Ec{g), defined as the minimum number of singlets (maximally 
entangled antisymmetric two-qubit states) which is needed to prepare the state g 
through local operations and classical communication |87| . In formulae: 

Ec{g) = lim ' . (27) 

n — >oo Tl 

The optimal convex decomposition of Eq. (|26|l can be found for symmetric states 
and turns out to be Gaussian, allowing for the determination of the entanglement of 
formation Ep'- 

Ep = max[0,/i(z>_)] , (28) 



with 



, , , {l+x)\ ( {l + xf \ {l~x)\ 
K^) = log2 —r— J— log2 



Ax \ Ax I Ax V Ax 



Such a quantity is, again, a decreasing function of , thus providing a quantification of 
the entanglement of symmetric states equivalent to the one provided by the logarithmic 
negativity £V- 

In the nonsymmctric case, an important result is that for any entangled two-mode 
Gaussian state the symplectic eigenvalue (and, consequently, the logarithmic 
negativity) can be estimated with remarkable accuracy by only determining the global 
purity Tr g^ of the state, and the two local purities Tr g\ 2 of each of the two reduced 
single-mode states gi — Tr^gi. For the aforementioned class of nonsymmctric thermal 
squeezed states, the estimate becomes actually an exact quantification, since the 
logarithmic negativity of these states is a function of the three purities only [3H| ■ These 
states are indeed the maximally entangled two-mode Gaussian states at fixed global 
and local purities, and thus they are the states one would like to produce and exploit 
in any continuous-variable quantum information processing. On the other hand, the 
symmetric instance, which carries the highest possible entanglement [221 (and which 
is the experimental product of the present paper) , is the one that enables continuous- 
variable teleportation of an unknown coherent state 1111 IT^ with a fidelity 
arbitrarily close to 1 even in the presence of noise (mixedness), provided that the 
state is squeezed enough (ideally, a unit fidelity requires infinite squeezing) . Actually, 
the fidelity of such an experiment, if the squeezed thermal states employed as shared 
resource are optimally produced, turns out to be itself a measure of entanglement and 
provides a direct, operative quantification of the entanglement of formation present in 
the resource jlH] . 



4. Effects of mode coupling on the entanglement generation 



As mentioned in the introduction, entanglement is very often produced by mixing on 
a beam splitter two squeezed modes. In the general case, the squeezed quadratures 



Entanglement of two-mode Gaussian states: characterization and experimental 



10 



have an arbitrary phase difference. We denote 6 + Tr/2 the phase difference between 
the two squeezed quadratures. The CM of the squeezed modes is then 



fa 

l/a 



0^ 




5c 
V c b' J 

while the CM of the two modes after the beam-splitter is 



where 



a sin^ 6 , b' = a cos^ 



sm^e 



k = 



cos^ 6* + (sin^ 6* + 1) 



sm 



2a 



I -a' 
2a 



k' 



n2 
a^ 



2a 



2a 



(29) 



ni 


k' 


k 


k' 


k' 


n2 


k' 


-k 


k 


k' 


ni 


k' 


k' 


~k 


k' 


n2 , 



(30) 



a I sm W cos f 



Let us first note that expression 1301 can be brought back to the expression given on 
fig. 5 of via local unitary operations which do not modify the entanglement. 

The CM of the squeezed (y4.± ) modes gives a good insight into the properties of 
the two-mode state. One can see that the intermodal blocks are zero, meaning that 
the two modes are uncorrelated. Consequently, they are the two most squeezed modes 
of the system (no further passive operation can select more squeezed quadratures). 
But one can also note that the two diagonal blocks are not diagonal simultaneously. 
This corresponds to the tilt angle of the squeezed quadrature. In order to maximize 
the entanglement, the two squeezed quadratures have to be made orthogonal, which 
can be done by a phase-shift of one mode relative to the other. 

It is easy in fact to compute the logarithmic negativity quantifying entanglement 
between the entangled modes A\ and Ai, when the two squeezed quadratures are 
rotated of 7r/2 + B. One has EjsjiJ- A^A^i) = -(1/2) logz>2, with 




cos2(6l) a4 + 6a2 + (a2 - 1)2 008(26*) + 1 L (31) 



The symplectic eigenvalue v is obviously a periodic function of 0, and it is globally 
minimized for = kir, with k £ Wi. The entanglement, in other words, is maximized 
for orthogonal modes in phase space, as already predicted in Rcf. [2]. Notice that 
this results holds for general nonsymmetric states, i.e. also in the case when the two 
modes Ai and A2 possess different individual squeezings. For symmetric states, the 
logarithmic negativity is depicted as a function of the single-mode squeezing a and 
the tilt angle 9 in figure ^ 

In the experiment we will discuss below, the entanglement is produced by a single 
device, a type-II OPO operated below threshold. When no coupling is present in the 



Entanglement of two-mode Gaussian states: characterization and experimental ... 11 




Figure 1. Logarithmic negativity as a function of tlie single-mode squeezing a 
and tlic tilt angle 8 between the two non-orthogonal quadratures in presence of 
mode coupling. 

optical cavity, the entangled modes are along the neutral axis of the crystal while 
the squeezed modes corresponds to the ±45° linear polarization basis. However, we 
have shown theoretically and experimentally in |2Uj that a coupling can be added via 
a birefringent plate which modifies the quantum properties of this device: the most 
squeezed quadratures are non-orthogonal with an angle depending on the plate angle. 
When the plate angle increases, the correlated (j4_) quadrature rotates of a tilt angle 
9 and the correlations are degraded. The evolution is depicted on figure [3 through the 
noise ellipse of the superposition modes. 




Figure 2. Frcsncl representation of the noise ellipse of the ±45° rotated modes 
when the coupling is increased. The noise ellipse of the —45° mode rotates and the 
noise reduction is degraded when the coupling increases while the +45° rotated 
mode is not affected. 

Eqn. 1311 shows that when coupling is present, it is necessary to perform an 
operation on the two modes in order to optimize the available entanglement. Before 
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Figure 3. Experimental setup. A continuous-wave frequency-doubled NdiYAG 
laser pumps below threshold a type II OPO with a A/4 plate inserted inside 
the cavity. The generated two-mode vacuum state is characterized by two 
simultaneous homodyne detections. The infrared output of the laser is used 
as local oscillator after filtering by a high-finesse cavity. The two couples 
{A/4, A/2} on each path are used to make arbitrary phase shift between orthogonal 
components of polarization. PD Lock: FND-100 photodiode for locking of the 
OPO. PD Split: split two-element InGaAs photodiode for tilt-locking of the 
filtering cavity. 

developing experimental measures of entanglement and optimization of the available 
resource in our system, let us detail our experimental setup. 

5. Experimental setup and homodyne measurement 

The experimental scheme is depicted on figure 13 and relies on a frequency-degenerate 
type-ll OPO below threshold. The system is equivalent to the one of the seminal 
experiment by Ou et al. but a A/4 birefringent plate has been inserted inside the 
optical cavity. When this plate is rotated, it results in a linear coupling between the 
signal and idler modes which induces above threshold a phase locking effect at exact 
frequency degeneracy |47[ I48| . This triply-resonant OPO is pumped below threshold 
with a continuous frequency-doubled Nd:YAG laser. The input flat mirror is directly 
coated on one face of the lOmm-long KTP crystal. The reflectivities for the input 
coupler are 95% for the pump (532nm) and almost 100% for the signal and idler 
beams (1064nm). The output coupler (R=38mm) is highly reflecting for the pump 
and its transmission is 5% for the infrared. At exact triple resonance, the oscillation 
threshold is less than 20 mW. The OPO is actively locked on the pump resonance 
by the Pound-Drever-Hall technique. The triple resonance is reached by adjustment 
of both the crystal temperature and the frequency of the pump laser. Under these 
conditions, the OPO can operate stably during more than one hour without mode- 
hopping. The birefringent plate inserted inside the cavity is exactly A/4 at 1064 nm 
and almost A at the 532 nm pump wavelength. Very small rotations of this plate 
around the cavity axis can be operated thanks to a piezoelectric actuator. 

Measurements of the quantum properties of arbitrary quadratures of light mode 
are generally made using homodyne detection 01]. When an intense local oscillator 
is used, one obtains a photocurrent which is proportional to the quantum noise of 
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the light in a quadrature defined by the phase-shift between the local oscillator and 
the beam measured. This photocurrent can be either sent to a spectrum analyzer 
which calculates the noise power spectrum or numerized for further treatments like 
tomographic measurements of the Wigner function 031 or selection 01]. As mentioned 
above, one can also characterize the entanglement by looking at linear combinations of 
the optical modes as opposed to linear combinations of the photocurrents |17lll8j . The 
two-modes which form the entangled state must be transformed via the beam splitter 
transformation, that is they are mixed on a 50/50 beam splitter or a polarizing beam- 
splitter into two modes which will be both squeezed if the original state is entangled. 

Homodyne detection allows for a simple and direct measurement of the 2x2 
diagonal blocks of the 4x4 CM. In order to measure the 2x2 off-diagonal 
blocks, linear combinations of the photocurrents can be used. One can thus obtain 
a complete measurement of the two-mode CM. However, in order to characterize 
simultaneously two modes a single phase reference is needed. To implement this, we 
have built a simultaneous double homodyne detection (Fig. |2| in box). The difference 
photocurrents are sent into two spectrum analyzers triggered by the same signal. 
Two birefringent plates inserted in the local oscillator path are rotated in order to 
compensate residual birefringence. A A/4 plate can be added on the beam exiting the 
OPO in order to transform the in-phase detections into in-quadrature ones (making 
the transformation {xi,pi,X2'P2) (•'£i,Pi,P2X2)). In such a configuration, two states 
of light with squeezing on orthogonal quadratures give in phase squeezing curves on 
the spectrum analyzers. 

6. Experimental measures of entanglement by the negativity 

As all experimental measurements, the measurement of the CM is subject to noise. 
It is thus critical to evaluate the influence of this noise on the entanglement. A 
quantitative analysis, relating the errors on the measured CM entries (in the A± 
basis) to the resulting error in the determination of the logarithmic negativity (the 
latter quantifying entanglement between the corresponding Ai and A2 modes) has 
been carried out and is summarized in Fig. ^ in absence of mode coupling. In general, 
the determination of the logarithmic negativity is much more sensitive to the errors 
on the diagonal 2x2 blocks a and /3 (referring to the reduced states of each mode, see 
Eq. dJ) of the CM r than on the off-diagonal ones (7, and its transpose 7"^, whose 
expectations are assumed to be null) . Let us remark that the relative stability of the 
logarithmic negativity with respect to the uncertainties on the off-diagonal block adds 
to the reliability of our experimental estimates of the entanglement. Notice also that, 
concerning the diagonal blocks, the errors on diagonal (standard form) entries turn out 
to affect the precision of the logarithmic negativity more than the errors on off-diagonal 
(non standard form) entries. This behavior is reversed in the off-diagonal block, for 
which errors on the off-diagonal (non standard form) entries affect the uncertainty on 
the entanglement more than errors on the diagonal (standard form) entries. 

Experimentally, we have measured the noise on the CM elements to be at best 
on the order of a few percents of the measured values for the diagonal blocks, 
corresponding to a fraction of a dB 021 ■ This is the case for the diagonal blocks 
which are well-known since they are directly related to the noise measurements of A^ 
and A-. The situation is less favorable for the off-diagonal blocks: the off-diagonal 
elements of these blocks show a large experimental noise leading which, as shown on 
Fig. E(b), may lead in some cases to unphysical CM, yielding for instance a negative 




(a) (b) 

Figure 4. Error SEj^f on the logarithmic negativity between modes Ai and A2, 
as a function of the error &T on the entries of the diagonal (a) and off-diagonal (b) 
2x2 blocks of the measured CM V in the A± basis, given by Eq. G^J. In plot (a): 
the solid red curve refers to equal errors (of value 5V) on the eight entries of the 
diagonal blocks, the dotted blue curve refers to equal errors on the four diagonal 
entries of the diagonal blocks while the dashed green curve refers to equal errors 
on the off-diagonal entries of the diagonal blocks (non standard form entries). At 
(5r > 0.16 some of the considered states get unphysical. In plot (b): the solid red 
curve refers to equal errors on the four entries of the off-diagonal block, the dotted 
blue curve refers to equal errors on the two off-diagonal entries of the off-diagonal 
block (non standard form entries), while the dashed green curve refers to equal 
errors on the diagonal entries of the off-diagonal block (standard form entries). 
At 5r > 0.19 some of the considered states get unphysical. 



determinant and complex values for the logarithmic negativity. In the following, we 
will set these terms to zero in agreement with the form of the CM of Eq. 123 

Let us first give an example of entanglement determination from measurements 
of CM elements, in the absence of mode coupling. Without the plate, the squeezing 
of the two superposition modes is expected on orthogonal quadratures (the ideal CM 
is then in the form I29| with 9 = 0). Spectrum analyzer traces while scanning the local 
oscillator phase are shown on figure [31 the rotated modes are squeezed on orthogonal 
quadratures. The state is produced directly in the standard form and the CM in the 
A± basis can be extracted from this measurement: 
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The resulting smallest symplectic eigenvalue is the geometric average of the two 
minimal diagonal elements : v- = 0.33, yielding a logarithmic negativity = 
— log2(£'-) = 1.60 between the modes Ai and A2. 



7. Experimental non standard form and optimization by linear optics 

As discussed previously, when the plate angle is increased, the state produced is not 
anymore in the standard form but rather similar to eqn. 1291 Figure El gives the 
normalized noise variances at 3.5 MHz of the rotated modes while scanning the local 
oscillator phase for an angle of the plate of 0.3°. The first plot shows that the squeezing 
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Figure 5. Normalized noise variances at 3.5 MHz of the ±45° modes while 
scanning the local oscillator phase. The first plot corresponds to in-phase 
homodyne detections and the second one in-quadraturc. Squeezing is well 
observed on orthogonal quadratures. (RBW fOO kHz, VBW 1 kHz) 



is not obtained on orthogonal quadratures. The CM takes the following form: 



r(p = o.3°) = 
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(0) (0) 9.54 -5.28 
V (0) (0) -5.28 3.45 
In this instance, the logarithmic negativity between Ai and A2 is much lower than 
the previous value: £V = 1.13. We will now study the action of passive operations 
on this two-mode state. 

As stated above, entanglement can be increased via passive operations performed 
simultaneously on the two modes. Such operations include phase-shifts and beam 
splitters, which can be readily performed on co-propagating, orthogonally polarized 
beams The minimal combination of waveplates can be shown to consist in three 
waveplates : two A/4 waveplates denoted Q and one A/2 waveplate denoted H. When 
using any combination of these three plates, the state can be put back into standard 
form which will maximize the entanglement. This operation consists in a phase-shift 
of the rotated modes. Figure El gives the normalized noise variances before and after 
this operation. The CM is changed to: 
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giving a logarithmic negativity £V = 1-32 between Ai and A2, larger than the value 
before the operation. It is also the maximal value than can be obtained considering 
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Figure 6. Normalized noise varianees at 3.5 MHz of the rotated modes wtiile 
scanning the loeal oscillator phase for an angle of the plate of 0.3°, before and 
after the non-local operation. The homodyne detections are in-quadrature. After 
this operation, squeezing is observed on orthogonal quadratures. 



the available entanglement. 

Let us remark again that this transformation is non-local in the sense of the EPR 
argument: it has to be performed before spatially separating the entangled modes for 
a quantum communication protocol for instance. 



8. Conclusion 

We have described the powerful tools underlying the description of continuous 
variables systems in quantum optics. These tools allow for a nice pictorial view 
of two-mode Gaussian entangled states. Specifically, we have illustrated their 
properties through the description and manipulation of the entanglement produced 
experimentally with an original device, a type-Il optical parametric oscillator 
containing a birefringent plate. We have demonstrated the capabilities of this system 
through entanglement measurements and manipulation of entangled states, showing, 
in particular, how entanglement can be maximized using purely passive operations. 

We have also studied quantitatively the influence of the noise, affecting the 
measurement of the elements of the CM, on the entanglement, showing that the most 
significant covariances (exhibiting the highest stability against noise) for an accurate 
entanglement quantification arc the diagonal terms of the diagonal single-mode blocks, 
and the off-diagonal terms of the intermodal off-diagonal block, the latter being the 
most difficult to measure with high precision. Alternative methods have been devised 
to tackle this problem |49l OHI based on direct measurements of global and local 
invariants of the CM. Such techniques have been implemented in the case of pulsed 
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beams |5(J| but no experiment to date has been performed for continuous-wave beams. 
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